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ABSTRACT 


An  experimental  study  was  performed  to  determine  the  feasibility  of 
a  laboratory  experiment  to  test  an  existing  theoretical  model  describing 
sound  propagation  into  a  fast  bottom  underlying  a  wedge  shaped  medium. 
Sand  under  fresh  water  was  found  to  satisfy  the  constraints  of  the  theor¬ 
etical  model  and  to  simulate  the  continental  shelf.  In  a  laboratory  ex¬ 
periment,  accuracy  of  density,  sound  speeds,  and  attenuation  was  shown 
to  be  sufficient  to  allow  quantitative  comparison  to  the  predictions  of 
the  beam  angle. 
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I.  INTRODUCTION 


The  behavior  of  radiation  in  a  wedge-shaped  medium  overlying  a  fast- 
bottom  medium  has  been  investigated  both  theoretically  and  experimentally 
with  both  optical  and  acoustic  sources.  For  example,  in  1971  Tien  and 
Martin  [Ref.  1]  examined  the  behavior  of  a  laser  beam  coupled  into  a  thin, 
tapered,  dielectric  film  deposited  on  a  substrate  with  a  higher  refrac¬ 
tive  index.  As  the  light  propagated  toward  the  apex  of  the  wedge  shaped 
film,  perfect  reflection  was  observed  until  the  changing  angle  of  inci¬ 
dence  decreased  below  the  critical  angle,  then  the  light  was  converted 
into  radiation  in  the  substrate. 

An  acoustic  analog  to  the  above  optical  study  was  performed  by 
Kuznetzov  [Ref.'  2]  in  1973.  Kuznetzov  developed  a  normal  mode  theory 
for  sound  propagation  in  both  a  wedge  and  underlying  half-space  substrate, 
where  the  substrate  sound  speed  was  faster  than  that  of  the  wedge. 
Kuznetzov  concluded  that: 

1.  Sound  traveling  toward  the  wedge  apex  would  be  totall  reflected 
until  the  angle  of  incidence  decreased  to  the  limiting  angle  of  total 
reflection. 

2.  Any  sound  incident  at  less  than  the  limiting  angle  of  total  re¬ 
flection  would  be  completely  refracted  into  the  underlying  half-space. 

This  total  refraction  would  occur  along  the  wedge/half-space  boundary 
from  the  wedge  apex  to  the  point  where  the  limiting  angle  was  first 
achieved. 

3.  Acoustic  energy  in  the  half-space  would  be  columnated  into  a 
well-defined  beam.  The  beam's  maximum  pressure  would  occur  at  an  angle 
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of  depression  (measured  from  the  plane  of  the  wedge/half-space  interface) 
that  would  lie  between  8  and  23,  where  6  is  the  wedge  angle. 

Kuznetzov  performed  a  series  of  experiments  that  supported  his  theory. 

The  formation  of  a  well-defined  beam  within  the  substrate  (for  the 
optical  case)  was  subsequently  demonstrated  by  Tien,  Smol insky  and 
Martin  [Ref.  3].  Reference  3  also  presented  two  theoretical  results: 

(1)  ray-optics  predicted  refraction  into  the  substrate  beyond  the  cutoff 
distance  (the  point  at  which  the  critical  angle  was  reached)  and  (2) 
theory  predicted  refraction  into  the  substrate  before  the  cutoff  distance. 

A  similar  well-defined  beam  was  predicted  and  observed  by  Si  gel  man, 
et  al.  [Ref.  4j  in  a  water-aluminum  system.  Maximum  pressure  occurred 
at  11  degrees  below  the  water-aluminum  interface  when  the  wedge  angle  was 
1.3  degrees.  This  was  well  outside  the  range  predicted  by  Kuznetzov. 

The  observed  beam  was  broader  than  predicted. 

In  1980,  Bradshaw  [Ref.  5]  extended  a  computer  model  first  developed 
by  Kawamura  and  Ioannou  [Ref.  6j  to  describe  the  behavior  of  sound  in  a 
fast  bottom  underlying  a  wedge  shaped  medium.  The  model  predicted  the 
formation  of  well-defined  beams  in  the  bottom,  and  predicted  the  effect 
of  attenuation  in  the  bottom  on  the  beam.  A  comparison  of  the  predic¬ 
tions  of  this  model  with  the  experimental  data  of  Netzorg  [Ref.  7]  gave 
qualitative  agreement,  but  since  the  experiment  did  not  match  all  model 
constraints,  quantitative  comparison  was  impossible.  It  is  the  purpose 
of  this  laboratory  analysis  to  determine  if  a  laboratory  experiment 
that  fulfills  the  constraints  of  the  model  can  be  designed.  If  at  all 
possible,  it  is  further  desired  to  model  the  continental  shelf  in  the 
design  of  the  experiment. 
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II.  CONSTRAINTS 


The  general  design  of  the  proposed  experiment  is  illustrated  in 
Figure  1.  Constraints  will  be  imposed  by  both  the  theoretical  model 
and  the  measurements. 

A.  MODEL  CONSTRAINTS 

The  Bradshaw  model  [Ref.  5]  requires  that  the  sound  speed  in  the 
bottom  exceed  that  in  the  overlying  wedge.  Physical  parameters  required 
to  be  known  are  the  density  and  sound  speed  of  the  two  media,  the  slope 
of  the  wedge,  and  the  attenuation  of  the  substrate  medium.  A  further 
constraint,  dependent  on  the  substrate  sound  speed,  will  be  discussed  in 
detail  later.  This  model  assumes  planar  waves  incident  on  the  wedge- 
bottom  interface. 

B.  OTHER  CONSTRAINTS 

Since  it  was  desired  to  model  a  real-world  environment,  specifically 
the  continental  shelf,  a  laboratory  set-up  consisting  of  fresh  water 
over  sand  was  an  obvious  first  choice.  The  sound  speed  of  fresh  water 
was  expected  to  be  about  1481  m/s  [Ref.  8]  and  the  speed  in  sand  about 
1700  m/s  [Ref.  9].  Thus,  the  requirement  of  a  faster  bottom  medium 
would  be  readily  fulfilled. 

A  wedge  of  water  would  require  a  high  frequency  (i.e.,  greater  than 
100  kHz)  to  permit  planar  incident  waves  in  a  laboratory-sized  experi¬ 
ment.  If  one  hundred  wavelengths  are  desired  to  ensure  plane  waves,  at 
100  kHz  the  source  would  need  to  be  1.5  m  away  from  the  wedge  apex. 


The  use  of  a  water/sand  system  would  permit  the  hydrophone  to  be 
moved  about  in  both  media;  the  use  of  plaster  to  simulate  rock,  for 
example,  would  not.  While  it  would  also  be  necessary  to  maintain  a  flat 
slope  of  at  least  five  degrees  between  the  two  media,  a  water-sand  sys¬ 
tem  was  expected  to  be  adequately  stable.  The  smoothness  of  the  slope 
is  required  by  the  model;  the  five  degree  slope  encompasses  the  maximum 
for  the  continental  shelf. 
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III.  EXPERIMENTAL  DESIGN 


A.  MATERIAL  SELECTED 

Fresh  (tap)  water  and  #30  fine  sand  were  the  media  used  in  the  experiment. 
The  grain  size  of  #30  fine  sand  varies  from  0.70  mm  to  0.15  mm.  It  was 
readily  demonstrated  that  the  sand  could  maintain  a  slope  as  large  as  25 
degrees.  To  remove  all  air  bubbles,  the  water/sand  system  was  heated 
with  immersion  heaters.  Bleach  was  added  at  a  ratio  of  one-half  gallon 
bleach  to  70  gallons  of  water  plus  sand  in  order  to  control  biologic 
growth. 

B.  VELOCITY  MEASUREMENT  TECHNIQUES 

An  excellent  summary  of  methods  for  measuring  sound  speeds  in  rocks 
and  minerals  is  given  by  Anderson  and  Liebermann  [Ref.  10]  .  Specialized 
methods  for  use  on  small  specimens,  iron  compounds,  etc.  or  for  high 
pressure  and  temperature  systems  were  eliminated  as  inappropriate  for 
the  experiment  being  considered.  Ultrasonic  interferometry  was  also  dis¬ 
carded  despite  its  exceptional  accuracy  (0.01  to  0.04  percent)  as  it  is 
best  applied  to  nonporous  samples  and  utilizes  frequencies  (10-60  MHz) 
that  would  cause  unacceptable  attenuation  in  the  verification  experi¬ 
ments.  Resonance  methods  with  appropriate  frequencies  and  accuracies 
were  undesirable  as  they  are  applicable  to  crystalline  specimens  only. 
Ultimately,  pulse  transmissions  were  selected  as  they  best  fit  the 
necessary  criteria.  This  method  is  applicable  to  either  fluid  or  por¬ 
ous  materials,  acceptable  over  a  frequency  range  of  50  kHz  to  10  MHz, 
and  accurate  to  within  one  to  three  percent. 


1.  Size  Criteria 


Steel -bound  glass  tanks  measuring  70  cm  square  by  60  cm  high 
were  readily  available.  Length  vs  depth  criteria,  from  Reference  10, 
require  that  the  layer  of  sand  on  the  tank  bottom  should  not  be  less 
than  14  cm  to  avoid  boundary  reflection  problems.  To  satisfy  this  con¬ 
straint,  a  sand  layer  of  20  cm  was  chosen. 

To  avoid  pulse  dispersion,  Reference  10  also  suggested  that  the 
sand  depth  not  be  less  than  five  wavelengths.  Since  the  sound  speed  of 
the  sand  was  expected  to  be  1700  m/s,  and  the  frequency  would  be  of  the 
order  of  100  kHz,  five  wavelengths  would  be  8.5  cm.  A  sand  depth  of 
20  cm  would  also  satisfy  this  restraint. 

To  avoid  particulate  scattering,  the  wavelength  must  be  at 
least  three  times  the  grain  size  [Ref.  10].  Since  the  largest  grain 
size  in  number  30  sand  is  0.070  cm,  a  wavelength,  equal  to  1.7  cm, 
readily  qualifies. 

Since  sound  speed  is  not  a  function  of  frequency  [Ref.  9],  re¬ 
sults  from  this  laboratory  analysis  were  expected  to  be  applicable  to 
frequencies  used  in  the  experiments. 

2.  Electronic  Equipment 

A  schematic  of  the  equipment  configuration  is  shown  in  Figure  2. 
All  components  were  off-the-shelf  and  readily  available. 

Output  from  a  General  Radio  model  1310  oscillator  with  a  fre¬ 
quency  range  of  2  Hz  to  2  MHz  was  fed  simultaneously  into  a  frequency 
counter  and  a  tone  burst  generator.  The  counter,  a  Hewlett-Packard 
5233L,  would  read  +10  Hz  at  100  kHz.  The  General  Radio  Type  396-A  tone 
burst  generator  was  used  to  generate  either  8  or  16  cycle  pulses  The 
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output  passed  through  a  Hewlett-Packard  HP  467-A  power  amplifier  before  being 
fed  to  a  transducer.  The  pulsed  signal  was  also  fed  to  the  trigger  in¬ 
put  of  a  Tektronix  type  Rm-503  oscilloscope. 

Three  types  of  transducers  were  used  as  sources.  The  first  was 
a  homemade,  mylar,  broad-band  transducer  with  a  3  cm  by  9  cm  active 
face.  This  source  was  quite  directional  at  100  kHz  but  proved  to  be 
difficult  to  move  through  the  sand  (for  sound  speed  and  attenuation 
measurements)  and  was  not  sufficiently  powerful  to  generate  usable  sig¬ 
nals  in  the  sand  (as  would  be  required  in  model  verification  runs).  The 
mylar  transducer  required  by  150  V  DC  supply  and  polarizing  network. 

A  second,  smaller,  homemade  transducer  was  tried  as  a  source. 

This  transducer,  a  matrix  of  ceramic  elements,  has  a  1  x  2  cm  active 
face  and  was  highly  directional  at  100  kHz.  This  transducer  offered 
little  improvement  over  the  first. 

Celesco  Industries  type  LC10  hydrophones  had  been  consistently 
used  as  receivers.  The  LC10  is  a  small  (0.97  cm  diameter  by  2.87  cm 
long)  cylinder,  with  a  receiving  range  of  0.1  to  120,000  Hertz.  The 
transducer  is  designed  to  be  omnidirectional  in  a  plane  perpendicular 
to  the  axis  of  the  cylinder  with  a  tolerance  of  +1  dB  at  100  kHz,  and 
omnidirectional  in  a  plane  containing  the  axis  of  the  cylinder  with  a 
tolerance  of  +2  dB  at  25  Hz.  The  omnidirectional  characteristic  was 
not  desirable  for  the  source,  as  this  would  increase  the  number  of  pos¬ 
sible  reflected  signals  from  the  walls,  water  surface,  etc.,  but  the 
small  size  and  higher  signal  strength  were  desirable.  Geometry  was  used 
to  isolate  desired  signals  and  will  be  described  when  appropriate. 
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Because  the  direct  and  reflected  pulses  arrived  from  different 
directions,  use  of  LClO's  as  the  source  and  receiver  for  attenuation 
measurements  initially  produced  non-reproducible  results.  Directivity 
measurements  in  the  plane  perpendicular  to  the  axis  o*  the  cylinder  were 
made  and  the  results  are  shown  in  Figs.  3  and  4.  Subsequent  measurements 
used  LC10  number  2338  as  a  source  and  LC10  number  2319  as  a  receiver. 

The  receiver  was  oriented  to  use  the  area  about  50  degrees  as  the  receiv¬ 
ing  face.  The  transducers  were  always  arranged  geometrically  to  use  the 
radial  plane. 

The  received  signals  were  amplified  20  dB  by  a  Hewlett-Packard 
HP-465A  amplifier,  then  passed  through  a  Spencer-Kennedy  Laboratories, 
Inc.  model  302  variable  electronic  filter  (set  at  60  kHz  high  pass)  to 
eliminate  low  frequency  mechanical  noise  present  in  the  laboratory  be¬ 
fore  being  passed  to  the  oscilloscope.  The  data  of  Fig.  5  demonstrates 
that  spherical  spreading  (1/r)  was  observed  for  all  LC10  source-receiver 
separations  greater  than  8.4  cm.  All  measurements  were  made  under  far- 
field  conditions. 
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Electronic  Equipment  Schematic 


Transmit 


Figure  3.  LC10  2338  Directivity  Pattern 
(units  of  volts) 
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IV.  MEASURED  DATA 


For  all  of  the  following  measurements,  the  sand  layer  in  the  tank 
was  kept  flat  and  horizontal,  as  only  physical  parameters  were  being  de¬ 
termined  and  not  wedge  characteristics. 

A.  DENSITY 

1 .  Water 

According  to  Lange's  "Handbook  of  Chemistry"  [Ref.  Ill,  the 
density  of  distilled  water  ranges  from  0.99913  g/cm3  at  15°C  to  0.99707 
g/cm3  at  25°C.  The  expected  density  of  room  temperature  water,  to  3 
significant  figures,  was  therefore  1.00  g/cm  .  Measured  volumes  of  the 
bleached  water  were  periodically  weighed,  and  the  observed  density  was 

3 

1.00  g/cm  ,  with  no  noted  variation. 

2.  Sand 

The  density  of  water  saturated  sand  was  measured  by  partially 
filling  a  weighed  100  ml  graduated  cylinder  with  a  water  over  sand  mix¬ 
ture,  observing  the  volume  of  each,  and  using  a  density  of  1.00  g/cm3 
for  the  water  to  calculate  the  weight  of  the  water  volume.  The  density 
of  the  water  saturated  sand,  for  eight  separate  measurements,  was 
1.98+0.03  g/cm3. 

B.  SOUND  SPEED 
1 .  Water 

Three  separate  sets  of  sound  speed  in  water  measurements  were 
made.  The  first  involved  the  use  of  the  mylar  transducer  as  the  source. 
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One  LCIO  was  placed  20.1  cm  away,  and  the  second  LC10  was  moved  along 
the  straight  line  defined  by  the  source  and  first  receiver.  The  time 
of  flight  between  the  two  receivers  was  measured.  A  simultaneous  check 
for  1/r  dependence  confirmed  these  measurements  were  all  in  the  far- 
field.  The  data  are  summarized  in  Table  1  of  Appendix  A.  The  average 
sound  speed  was  1454  m/s  with  an  (N-l)  variance  of  33  m/s. 

The  second  set  of  data  utilized  a  reflected  signal  to  increase 
the  total  path  length.  The  mylar  source  was  fixed,  and  the  receiver, 
placed  between  the  source  and  the  wall,  was  moved  along  a  perpendicular 
line  between  the  source  and  the  wall.  The  time  of  flight  was  the  time 
difference  between  the  arrival  of  the  direct  and  reflected  signals;  the 
distance  was  twice  the  distance  from  the  receiver  to  the  wall.  These  re¬ 
sults  are  summarized  in  Table  2  of  Appendix  A.  The  average  speed  was 
1469  m/s,  with  a  variance  of  14.7  m/s. 

In  the  third  set  of  measurements,  the  distance  over  which  the 
time  of  flight  was  measured  was  fixed.  The  LC10  source  was  suspended 
11.0  cm  from  the  tank  wall,  the  second  LC10,  placed  further  from  the 
wall  than  the  source,  was  moved  along  the  perpendicular  formed  by  the 
wall  and  the  source.  As  the  receiver  was  moved,  the  distance  differ¬ 
ence  between  the  direct  path  signal  and  the  signal  reflected  from  the 
wall  remained  a  constant  22.0  cm.  The  time  difference  between  the  ar¬ 
rival  of  the  two  signals  was  recorded.  Table  3  of  Appendix  A  summarizes 
the  results.  The  average  speed  was  1460  m/s,  the  variance  was  19  m/s. 

The  data  compared  well  with  theory.  From  Kinsler  and  Frey 
[Ref.  8],  the  equation  of  sound  for  distilled  water  is  given  by 

cw  =  1403  +  5t  -  0.06t2  +  0.0003t3  (1) 
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where  t  is  the  temperature  in  degrees  Celcius.  A  comparison  of  measured 
and  theoretical  values  is  presented  in  Table  VI-1.  Note  that  direct 
path  results  showed  greater  variation,  but  all  results  were  within  two 
percent  of  expected  values. 

2.  Sand 

Five  separate  sets  of  sound  speed  measurements  were  made  in  the 
sand;  four  in  conjunction  with  attenuation  data,  and  one  to  check  for  a 
sound  speed  gradient.  All  three  source  transducers  were  used.  For  the 
attenuation  runs,  the  source  was  buried  at  a  fixed  depth,  probing  was 
with  an  LC10  receiver  inside  a  glass  tube  for  the  depth  of  maximum  signal 
amplitude.  The  glass  tubes  were  necessary  to  prevent  abrasive  damage  to 
the  receiver.  The  data  are  summarized  in  Tables  1  through  4  of  Appendix 
8. 

The  fifth  set  of  data  utilized  the  two  LClO's  fixed  30  cm  apart. 
The  transducers  were  initially  placed  in  water  and  sand  was  added  to 
progressively  bury  them.  The  entire  process  was  carried  out  underwater, 
so  degassed,  saturated  sand  was  used.  The  data  are  summarized  in  Table 
5  of  Appendix  B.  Table  IV-2  lists  the  resultant  sound  speeds. 

a.  Theoretical  Sound  Speed  in  Sand 

Since  a  search  of  the  literature  revealed  many  theoretical 
and  experimental  values  for  the  speed  of  sound  in  sand  under  brine,  but 
not  under  fresh  water,  and  since  the  values  given  were  greater  than 
those  observed  in  this  experiment,  the  following  analysis  was  performed 
to  derive  a  theoretical  value  for  sound  speed  in  sand  under  fresh  water. 

From  Urick  [Ref.  12],  the  density  of  a  liquid-sediment 
mixture  is  given  by 
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Table  IV-1.  Measured  vs  Theoretical  Sound  Speed  in  Water 
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p.  =YP  +  (1  -  v)p 

Mrmx  '  w  v  "  s 

where  pmi.  =  the  density  of  the  mixture, 


Y  =  the  porosity, 

p. ,  =  density  of  water. 


and 


p$  =  the  density  of  individual  sand  grains. 


From  equation  (2),  using  the  measured  value  pmix  =  1.98  g/cnr  and 

3 

p$  =  2.69  g/cm  from  Reference  9,  gives  y  =  0.42.  Also  from  Reference 
12,  the  speed  of  sound  in  the  saturated  sediment  is  given  by 


cmix  =  ftv  +  ps(1  *  Y))(kwY  +  M1  -  y))J 


-1/2 


(3) 


where  k  and  k  are  the  compressibility  of  the  water  and  sand,  respec- 
w  s 

tively.  Values  for  kw  and  k$  were  calculated  from  the  relationships 


cw  =  W1'2 

(4) 

and 

s  - 

(5) 

Setting  cw  =  1481  m/s  [Ref.  8],  and  pw  =  1.0  g/cm'3,  and  solving  for 
kw  yields 

kw  =  4.56  x  IQ'13  m  sec2  g'1  (6) 

The  sand  in  the  experiment  was  composed  of  quartz  particles.  From 
Anderson  and  Lieberman  [Ref.  10],  page  360,  c$  was  set  to  7000  m/s. 
Solving  equation  (5)  yields 

ks  =  7.59  x  IQ'15  m  sec2  g"1  (7) 

Substituting  the  above  values  in  equation  (3)  and  solving  for  cm^x 

yields  an  expected  speed  for  sand  under  fresh  water  of  1605  m/s. 

A  comparison  of  the  expected  value  with  the  sound  speeds 
summarized  in  Table  IV-2  shows  close  correlation  with  the  value  obtained 


X 
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TABLE  IV-2.  Measured  Sound  Speeds  in  Sand 


Measurement 

Method 

Sound  Speed 
(m/s) 

(N-l)  a 
(m/s) 

Variable  distance 

1579 

58.6 

Variable  distance 

1597 

44.7 

Variable  distance 

1434 

71.9 

Variable  distance 

1567 

32.6 

Fixed  source/receiver 

1613 

0.0 
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using  a  fixed  source-receiver  distance  (1613  m/s)  but  not  close  agreement 
with  the  majority  of  values  obtained  from  measurements  utilizing  glass 
tubes.  The  use  of  glass  tubes  disturbs  the  medium,  and  the  distance  from 
source  to  receiver  is  not  accurately  measurable.  Over  the  short  distances 
used,  a  0.5  cm  error  in  the  distance  would  lead  to  significant  errors. 

Thus,  the  value  obtained  from  the  fixed  source  and  receiver  method  is  more 
likely  accurate,  and,  not  unexpectedly,  shows  closer  agreement  with  the 
theoretical  values. 

b.  Sound  Speed  Gradient 

As  illustrated  in  Fig.  6,  Hamilton  [Ref.  13]  predicts  a  pro¬ 
nounced  sound  speed  gradient  in  a  fine  sand  sediment.  A  change  in  slope 
of  27  s"1  is  specified  at  zero  depth.  Although  the  gradients  were  estab¬ 
lished  for  sand  under  brine,  they  are  attributed  to  a  decrease  in  sediment 
porosity  with  pressure,  to  temperature  effects,  and  to  pressure  effects  on 
the  pore  water,  and  would  therefore  be  expected  to  be  present  in  fresh 
water  systems  as  well.  Sound  speed  measurements  tabulated  in  Table  5  of 
Appendix  B  were  taken  specifically  to  check  for  the  presence  of  a  gradient. 
None  was  found  within  the  accuracy  limitations  of  the  method.  For  a  sound 
speed  of  1613  m/s,  a  2 %  error  (see  Sec.  1 1 1 -A )  would  be  32  m/s  and  would 
mask  a  27  m/s  change.  Even  if  a  27  m/s  gradient  were  present,  the  re¬ 
sulting  62  m  radius  of  curvature  would  not  be  sufficient  to  be  detectable 
over  the  experimental  ranges.  The  lack  of  a  detectable  gradient,  although 
not  mirroring  actual  continental  shelf  conditions,  is  consistent  with  the 
isovelocity  substrate  assumed  in  Bradshaw's  model  [Ref.  5j. 
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OBSERVED  Cs  =  16 13  Z  0 


Figure  6.  Sound  Speed  Gradient  in  Brine 
Saturated  Sand  (adapted  from 
reference  13.) 


C.  ATTENUATION 

1 .  Model  Limitations 

It  was  mentioned  in  !><■:  cion  1 1 -A  that  a  further  constraint,  to  be 
defined  once  the  substrate  sound  speed  had  been  determined,  must  be  ful¬ 
filled  for  verification  of  Bradshaw's  model.  This  constraint  is  that  the 
signal  must  propagate  beyond  a  specific  range  determined  as  follows.  The 
model  uses  a  Green's  function  analysis  to  obtain  t. ;e  pressure  and  phase 
distribution  of  the  sound  in  the  bottom.  To  solve  the  appropriate  Green's 
function,  it  was  necessary  to  assume  that 

k2jr  -  r'|  »  2 it  (8) 

where  k2  =  u/cs  is  the  wave  number  for  the  sound  in  the  sand.  Using  the 
determined  sound  speed  of  1613  m/s  and  assuming  a  frequency  of  100  kHz, 
k^  equals  389  m”^.  Thus,  it  is  required  that  jr  -  r'|  »  0.016  m 

The  quantity  jr  -  r"|  is  defined  in  Fig.  7.  The  quantity  r'  is 
the  distance  from  the  wedge  apex  to  a  point  on  the  wedge-substrate 
boundary  (integration  in  the  model  is  carried  out  over  r').  The  quantity 
r  is  the  distance  from  a  point  along  the  wedge-substrate  boundary,  called 
X,  to  the  field  point.  The  distance  X  is  the  distance  at  which  the  cri¬ 
tical  angle  is  reached  and  where  lowest  mode  transmission  into  the  sub¬ 
strate  begins.  The  distance  is  defined  by 


w 

4  f  sin6c  sin6’ 

(9) 

=  arccos  (cw/cs) 

(10) 

where  f  is  the  frequency, 

3  is  the  wedge  angle 

cw,  c$  are  the  sound  speeds  in  the  wedge  fluid  and  substrate, 
respectively, 
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and  6C  is  the  critical  angle. 

For  f  =  100  kHz,  a  wedge  angle  of  3°,  c  *  1469  m/s  and  c  =  1613  m/s,  X 

W  9 

would  equal  0.17  m. 

Typically,  within  the  model,  the  integration  limit  is  from  nX  to 
NX,  where  n  =  0,1,2,...,  N  =  1,2,3,...,  N>n.  In  the  simplest  case,  the 
integration  would  be  carried  out  from  0  to  X.  This  places  a  constraint 
on  r,  and  thus  on  permitted  attenuation. 

Consider,  as  a  minimum,  a  factor  of  10  to  satisfy  the  "much 
greater  than"  inequality;  this  would  require  |r  -  r'\  to  be  greater  than 
0.16  m.  This  implies  mathematically  that  no  measurements  may  be  made  at 
point  r  within  a  circle  of  radius  16  cm  about  the  point  r",  and  since  r' 
can  vary  from  0.0  to  0.17  m,  an  area  of  "forbidden  territory"  is  traced 
about  the  wedge  apex  and  boundary  surface.  Thus,  to  test  the  model, 
since  measurements  must  be  made  outside  this  area,  the  attenuation  must 
not  be  so  great  that  insufficient  signal  remains  beyond  the  forbidden  zone. 
(If  "much  greater  than"  were  taken  to  be  100  times,  measurements  would 
have  to  be  made  over  1.6  m  from  the  boundary). 

2.  Results 

The  methodology  for  collecting  attenuation  data  was  described  in 
Section  IV-B-2.  The  data  were  processed  graphically.  Taking  the 
natural  logarithm  of  the  well-known  equation 

V  -a  r 

V  =  ?°e  (11) 

where  V  is  the  measured  voltage, 

VQ  is  the  source  voltage, 
r  is  the  distance  from  the  source, 
and  a  is  the  attenuation. 
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the  linear  relation 


ln(Vr)  =  ln(VQ)  -  ar  (12) 

is  obtained.  Thus  a,  in  nepers  per  meter,  will  be  the  slope  of  a  graph 
of  ln(Vr)  vs  r.  Graphs  of  the  four  data  sets  of  Appendix  B  are  shown  in 
Figures  8  through  11.  A  summary  of  the  attenuation  values  is  presented 
in  Table  IV-3. 

Figure  12  shows  a  plot  of  the  attenuation  vs  frequency.  The  line 
indicates  the  slope  of  any  line  having  a  dependence  of  attenuation  on  the 
first  power  of  frequency.  All  data  fall  within  the  limits  of  data  col¬ 
lected  by  Hamilton  [Ref.  14]. 

3.  Reflection  Coefficients  and  Signal  Level 

Normal  incidence  pressure  reflection  and  transmission  coefficients 
were  measured  as  an  accuracy  check  on  the  measured  p  and  c  values,  and 
also  to  provide  an  estimate  of  the  sound  pressure  level  that  would  be 
transmitted  from  the  wedge  to  the  underlying  medium.  A  summary  of  the 
measured  data  is  presented  in  Appendix  C. 

The  reflection  coefficient  R  is  calculated  from 


R  =  '  pwcw  ~  pscs ' 


!pwcw+  pscs! 


and  the  transmission  coefficient  T  from 


(13) 


T  =  1  -  R  (14) 

3  3 

With  ow  =  1.00  g/cm  ,  p$  =  1.98  g/cm  ,  cw  =  1469  m/s,  and  c$  *  1613  m/s, 
R  is  calculated  to  be  0.37  and  T  to  be  0.63.  The  measured  values  of  re¬ 
flection  were  0.31  and  0.29.  The  measured  transmission  coefficient  was 

0.65.  These  values  indicate  that  the  measured  density  and  sound  speed 
values  are  within  the  expected  range. 
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Figure  10.  Attenuation  in  Sand 


100 

Frequency  (kHz) 


igure  12.  Attenuation  vs  Frequency 


From  Table  1  of  Appendix  C,  a  received  voltage  of  6  volts  is 
typical.  Assuming  a  desired  reading  on  the  order  of  1  volt  for  test 
runs,  a  drop  in  level  of  no  more  than  15.6  dB  is  indicated.  At  100  kHz, 
the  measured  attenuation  of  56.6  dB  per  meter  would  indicate  measure¬ 
ments  could  be  made  at  least  27  cm  from  the  wedge  boundary.  This  satis¬ 
fies  the  limit  imposed  on  r  in  Section  IV-C-1. 


t 

\ 
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TA8LE  1V-3.  Attenuation  Results 


f (kHz) 

a(Ne/m) 

ct(dB/m) 

Correlation 

Coefficient 

75.10 

7.26 

63.2 

-0.90 

100.00 

6.50 

66.6 

-0.99 

100.04 

7.39 

64.3 

-0.76 

130.00 

9.79 

85.2 

-0.92 
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V.  ERROR  ANALYSIS 


It  is  important  to  know  if  the  physical  properties  were  determined 
with  sufficient  accuracy  to  verify  Bradshaw's  model.  Bradshaw  developed 
an  analytical  equation  for  the  depression  angle  of  the  beam  formed  in 
the  substrate.  As  this  analytical  equation  was  based  on  model  results, 
and  since  the  equation  is  dependent  on  all  measured  physical  properties, 
the  equation  should  provide  a  realistic  basis  for  error  analysis.  The 
analytical  equation  is 


0D  =  17.22 


5°-329  ec°-”2  (Vos)‘°-254 


05) 


where  9g  is  the  beam  depression  angle  in  degrees, 

3  is  the  wedge  angle  in  degrees, 

9C  is  the  critical  angle  in  radians, 
and  pw/ps  is  the  ratio  of  the  density  of  water  to  the  density  of 
saturated  sand. 

The  critical  angle  is  defined  by 

9C  =  arccos  (cw/c$)  (16) 

where  c,  and  c.  are  the  sound  speeds  in  water  and  sand,  respectively, 
w  s 

Using  standard  error  analysis,  and  noting  pw/p$  =  l/p  , 


deD«J'Vps)  ■  S“ dB  ♦  i°  d6c +  i' " dp 


(17) 


=  5.66  9C  0,772  (l/ps)  '  0,25480,671dB 
+  13.2980'329  (1/ps)  '  O,2549cO,228d0c 
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+  4.37s0’329  e°-772P°-74V 

c  s  s 

Similarly 

de„  =  50c  4c  *  56c  dc.  (18) 

c  ^  "  ^7  5 

»[^dcs-i_dcw^l  -  (c„/cs)2]-'« 

where  cw  =  1469*15  m/s  from  Table  2.  An  error  of  3  percent*  in  keeping 
with  the  worst  case  predicted  accuracy  from  Reference  10,  was  assigned  to 
sound  speed  in  saturated  sand.  Therefore,  cs  =  1613*48  m/s  and  d6c  *  0.04 
radians. 

Assuming  the  wedge  angle  can  be  measured  to  within  0.5  degrees,  the 
following  values  were  used  to  evaluate  equation  17. 

Ps  =  1.98  g/cm3  dp$  =  0.03  g/cm3 

3=3°  d8  =  0.5° 

0C  =  0.4257  rad  d0c  *  0.04  rad 

Therefore 

dSD  =  3.64  +  0.75  +  0.16 
=  4.55° 

o 

9D  for  the  above  parameters  is  15.21  degrees.  The  deviation  of  0.75 
attributable  to  sound  speed  measurements  is  roughly  5  percent,  but  is 
dwarfed  by  the  3.64  (24  percent)  error  attributable  to  the  wedge  angle 
measurement.  A  method  of  accurately  measuring  and  controlling  the 
wedge  angle  will  be  a  necessity  for  successful  model  test  runs. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  laboratory  sized  experiment  to  test  the  predictions  of  Bradshaw's 
model  is  possible,  using  off-the-shelf  equipment  in  a  water-over-sand 
system.  Measurements  can  be  made  of  all  necessary  physical  parameters, 
and  attenuation  is  not  so  severe  as  to  preclude  measurements  at  suffi¬ 
cient  distances  to  satisfy  model  constraints. 

Densities  and  sound  speeds  were  obtained  that  agreed  with  literature 

or  theoretical  values.  The  densities  of  water  and  saturated  sand  were 
3  3 

1.00  g/cm  and  1.98  g/cm  ,  respectively.  The  sound  speed  in  water  was 
1469  m/s;  the  sound  speed  in  saturated  sand  was  1613  m/s. 

The  use  of  a  rigid  array  of  receivers  is  highly  recommended  for  all 
measurements  within  the  sand  in  order  to  control  source-receiver  spacing. 

Model  results  are  strongly  dependent  on  the  wedge  angle.  The  sand 
will  maintain  a  flat,  sloped  surface,  but  great  care  in  the  formation  and 
measurement  of  the  wedge  angle  will  be  required. 

It  was  found  necessary  to  maintain  a  high  bleach  concentration  in 
order  to  control  gas  producing  biologic  growth.  Every  change  of  water 
should  be  accompanied  by  a  renewal  of  the  bleach  level. 
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APPENDIX  A 

SOUND  SPEED  IN  WATER  -  DATA 


TA8LE  1.  Two  Receivers  Direct  Path  Data 

Source  Transducer  =  Homemade  Mylar 


Separation 

(cm) 

Time 

O0~5  sec) 

Sound  Speed 
(meters/second) 

6.8 

4.50 

1511 

11.7 

7.95 

1472 

16.7 

11.70 

1427 

19.0 

13.30 

1428 

21.4 

14.90 

1436 

26.3 

18.10 

1453 

c  =  1454  meters/second 
a  *  32.6  [2.2%) 


TABLE  2.  Reflected  Path,  Varied  Distance  -  Data 
Source  Transducer  =  Homemade  Mylar 


Separation 

(cm) 

Time 

(10*^  sec) 

Sound  Speed 
(meters/second) 

28.8 

19.55 

1473 

38.8 

26.32 

1474 

48.8 

33.25 

1468 

58.8 

40.15 

1464 

68.8 

47.85 

1438 

78.8 

53.70 

1467 

88.8 

59.60 

1490 

98.8 

67.00 

1475 

c  *  1469  meters/second 
a  =  14.7  (1.0%) 


1 
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TABLE  3.  Reflected  Path,  Fixed  Distance  -  Data 


Source  Transducer  =  LC10  #2338 


Separation 

(cm) 

Time 

(10“®  sec) 

Sound  Speed 
(meters/second) 

22.0 

14.92 

1474 

22.0 

15.32 

1436 

22.0 

14.88 

1478 

22.0 

14.92 

1474 

22.0 

14.92 

1474 

22.0 

14.92 

1474 

22.0 

15.32 

1436 

22.0 

15.32 

1436 

22.0 

15.08 

1459 

c  =  1460  meters/second 
o  =  18.8  (1.3%) 


APPENDIX  8 

SOUNO  SPEED  AND  ATTENUATION  IN  SAND  -  DATA 
TABLE  1.  Sound  Speed  and  Attenuation  in  Sand 


Source  Transducer  =  Homemade  Mylar 

Frequency  =  100.04  kHz 

Source  placed  10  cm  below  sand  surface. 


Distance 

Time 

Amplitude 

Sound  Speed 

(cm) 

(10-5  sec) 

(mV) 

(meters/second) 

15.15 

9.04 

11.0 

1676 

16.15 

10.10 

19.0  (18.0)* 

1599 

19.15 

12.25 

7.0  (  7.2) 

1563 

19.60 

12.25 

9.0 

1600 

22.05 

13.80 

7.0  {  7.2) 

1598 

23.75 

15.40 

3.0  (  2.8) 

1542 

24.10 

17.55 

7.5 

1373** 

25.55 

15.95 

4.0  (  4.0) 

1602 

27.65 

17.00 

3.5  (  3.2) 

1626 

28.15 

19.15 

3.0  (  2.4) 

1470 

29.10 

18.05 

3.0 

1612 

29.70 

18.10 

6.0 

1641 

31.90 

20.20 

1.0 

1579 

31.90 

21.25 

1.5 

1501 

34.20 

22.85 

1.5 

1497 

38.00 

no  signal 

— 

— 

♦changed  scale 

c  =  1579  meters/second 

♦♦omitted  point 
(outside  2a) 

a  =  58.6  (3.7*) 

TABLE  2.  Sound  Speed  and  Attenuation  in  Sand 

Source  Transducer  =  Homemade  Mylar 
Frequency  =  100.00  kHz 

Source  placed  14.5  cm  below  surface  of  sand. 


Distance 

(cm) 

Time 

(10-5  sec) 

Amplitude 

(mV) 

Sound  Speed 
(meters/second) 

10.30 

6.38 

27.0 

1614 

18.39 

11.60 

10.0 

1585 

19.15 

12.12 

42.0* 

1580 

23.71 

14.26 

4.5 

1663 

29.22 

18.94 

2.8 

1543 

♦point  omitted  c  =  1597  meters/second 

for  attenuation 

a  *  44.7  (2.8%) 
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TABLE  3.  Sound  Speed  and  Attenuation  in  Sand 

Source  Transducer  =  Homemade  Ceramic 
Frequency  =75.1  kHz 

Source  placed  13.0  cm  below  sand  surface. 


Distance 

(cm) 

Time 

(10'5  sec) 

Amplitude 

(mV) 

Sound  Speed 
(meters/second) 

12.2 

9.05 

25.0 

1348 

12.5 

8.00 

20.0 

1562 

23.5 

16.50 

8.0 

1424 

20.6 

14.35 

3.0 

1436 

27 .4 

19.70 

2.5 

1391 

32.3 

22.35 

2.5 

1445 

c  =  1434  meters/second 
a  -  71.9  (5.0%) 
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TABLE  4.  Sound  Speed  and  Attenuation  in  Sand 


Source  Transducer  =  Homemade  Mylar 

Frequency  *  130.00  kHz 

Source  placed  7.5  cm  below  sand  surface. 


Distance 

(cm) 

Time 

(10'5  sec) 

Amplitude 

(mV) 

Sound  Speed 
(meters/second) 

12.0 

7.7 

10.0 

1558 

13.0 

8.0 

— 

1625 

18.7 

12.0 

3.0 

1565 

23.0 

14.9 

1.8 

1544 

26.1 

16.5 

0.8 

1582 

27.3 

18.1 

0.6 

1508 

31.5 

20.2 

0.6 

1559 

32.2 

20.2 

0.6 

1594 

32.6 

20.8 

0.4 

1571 

c  =  1567  meters/second 


TABLE  5.  Sound  Speed  Gradient  Data 


Source 

Transducer  = 

IC10 

#2338 

Frequency  =  120.05 

kHz 

Source 

and  receiver 

fixed 

at  30.0  cm  separation. 

Water 

(cm) 

Sand 

(cm) 

Time 

(IQ'5  sec) 

Amplitude 

(V) 

Sound  Speed 
(meters/second) 

11.0 

---- 

20.2 

5.4 

1485 

27.0 

— — 

20.2 

10.4 

1485 

30.0 

0.0 

18.6 

0.5 

1613 

29.0 

1.0 

18.6 

1.0 

1613 

27.0 

3.0 

18.6 

2.5 

1613 

r 

1 

25.5 

4.5 

18.6 

1.35 

1613 

24.0 

6.0 

18.6 

1.35 

1613 

23.0 

7.0 

18.6 

1.35 

1613 

21.0 

9.0 

18,6 

1.35 

1613 

‘ 

20.0 

10.0 

18.6 

1.35 

1613 

f 

17.5 

12.5 

18.6 

1.35 

1613 

15.0 

15.0 

18.6 

1.35 

1613 

12.5 

17.5 

18.6 

1.35 

1613 

«  '  1 
i: 

10.0 

20.0 

18.6 

1.35 

1613 

APPENDIX  C 


REFLECTION  AND  TRANSMISSION  COEFFICIENTS  -  DATA 

TABLE  1.  Surface  Reflection  Coefficient 


Direct 

Path 

Reflected 

Path 

Reflection  Coeff. 

Distance 

Amp. 

Distance 

Amp. 

(cm) 

(V) 

(cm) 

(V) 

22.5 

5.8 

48.5 

2.9 

1.077 

30.0 

5.4 

37.2 

4.0 

0.9185 

2  0.0 

6.2 

40.0 

3.3* 

1.010 

20.0 

6.4 

40.0 

3.0** 

0.9880 

R  = 

0.998 

a  = 

0.06 

*  Directivity  (.241/. 254) 

**  Directivity  (.254/. 241) 

*  and  **  measured  at  normal  incidence 
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TABLE  2.  Bottom  Reflection  Coefficient,  Normal  Incidence 


K 


k 

i 


Distance 

Direct  =  20.0  cm 
Reflected  =  40.0  cm 


Directivity  Index 
Direct  =  0.254 
Reflected  *  0.241 


Amplitude  (V) 

Reflection 

Coefficient 

Di rect 

Reflected 

5.4 

0.8 

0.31 

5.7 

0.8 

0.30 

6.5 

0.8 

0.26 

7.2 

1.0 

0.29 

7.1 

1.0 

0.30 

7.1 

0.9 

0.27 

7.0 

1.0 

0.30 

7.0 

0.9 

0.27 

7.2 

1.2 

0.35 

10.0 

1.2 

0.25 

8.6 

1.1 

0.27 

R  *  0.29 
a  =  0.028 
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TABLE  3.  Bottom  Reflection  Coefficient,  Normal  Incidence 


Distance 

Direct  =  19.2  cm 
Reflected  =  40.4  cm 


Directivity  Index 
Direct  =  0.254 
Reflected  =  0.241 


Amplitude  (V) 


Di rect 

Reflected 

7.2 

1.1 

7.0 

1.2 

7.5 

1.25 

7.0 

1.2 

7.1 

1.0 

7.6 

1.0 

7.6 

1.0 

7.5 

0.8 

6.5 

0.8 

5.6 

0.8 

6.6 

0.9 

7.4 

0.8 

7.4 

0.95 

7.5 

1.0 

7.0 

1.0 

7.0 

1.1 

7.5 

1.2 

7.0 

1.1 

Reflection 

Coefficient 

0.34 

0.38 

0.37 

0.38 

0.31 

0.29  R  =  0.31 

0.29  o  =  0.044 

0.23 
0.27 
0.32 
0.30 
0.24 
0.28 
0.30 
0.32 
0.35 
0.35 
0.35 
0.29 


7.0 


0.9 


53 


TABLE  4.  Transmission  Coefficient,  Normal  Incidence 


Direct  Path  =  20.0  cm  in  water 

Transmitted  =  10.0  cm  in  water,  10.0  cm  in  sand 


Frequency 

(kHz) 

Direct 

Amplitude 

(volts) 

Transmitted 

Transmission 

Coefficient 

125 

8.7 

6.8 

0.78 

130 

9.0 

7.0 

0.78 

125 

8.7 

5.0 

0.57 

120 

7.8 

5.0 

0.68 

115 

7.4 

3.6 

0.49 

no 

6.9 

3.6 

0.52 

120 

7.8 

5.4 

0.69 

125 

8.7 

5.2 

0.60 

130 

9.0 

7.6 

0.84 

125 

8.7 

5.2 

0.60 

120 

7.8 

5.2 

0.67 

115 

7.4 

3.8 

0.51 

115 

7.4 

3.8 

0.51 

120 

7.8 

5.2 

0.67 

125 

8.7 

5.3 

0.61 

130 

9.0 

7.6 

0.84 

T  =  0.65 

J 

j;  a  =  0.116 

i 

1 1 
r  " 
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